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^ | Abstract 

{Sj ■ We offer a new Hamiltonian formulation of the classical Pais-Uhlenbeck Oscillator and 

consider its canonical quantization. We show that for the non-degenerate case where the 
frequencies differ, the quantum Hamiltonian operator is a Hermitian operator with a positive 
spectrum, i.e., the quantum system is both stable and unitary. A consistent description of the 
D \ degenerate case based on a Hamiltonian that is quadratic in momenta requires its analytic 

continuation into a complex Hamiltonian system possessing a generalized "P7~-symmetry (an 
involutive antilinear symmetry). We devise a real description of this complex system, derive 
an integral of motion for it, and explore its quantization. 
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X ■ 1 Introduction 

a: 

Pais-Uhlenbeck (PU) Oscillator p] is the simplest and by far the best known toy model for a higher 
derivative theory. The interest in this model is motivated by the fact that it provides an opportunity 
to explore the possibility of solving a long standing problem related with the non-unitarity of higher 
derivative theories (of gravity). These theories are particularly interesting because they are known 
to be perturbatively renormalizable [2]. 

Consider the following equation of motion [3]. 

+az^ + f3z = 0, (1) 

where z is a real- valued function of time, z^ denotes the fc-th derivative of z, and a and (3 are real 
and positive parameters related to a pair of frequencies U\ and u 2 according to 

a := uj\ + lo%, (3 := cu^o;^. (2) 
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Throughout this article we define the classical PU oscillator as the classical dynamical system 
determined by the forth-order equation of motion (CQ). Our aim is to obtain a unitary and stable 
quantum system that has the classical system defined by ([1]) as its classical limit. Here by unitarity 
and stability of a quantum system we mean that the corresponding quantum Hamiltonian operator 
is Hermitian and its spectrum is bounded from below. 

It is well-known that (JTJ can be derived from a Lagrangian of the form [T| [3] 

L = ^f-az 2 + Pz 2 ), (3) 

where \x is an arbitrary positive real mass parameter and each over-dot stands for a time-derivative. 
Introducing a new coordinate variable x according to x — z, one can turn ([3]) into the Lagrangian of 
a constraint second derivative theory and apply the machinery of Dirac's constraint quantization to 
arrive at a Hamiltonian description of the PU Oscillator based on the following classical Hamiltonian 

H = ^ + x Pz + ^-^. (4 ) 
2fi yz 2 2 v ; 

In Ref. [5], the authors describe how the standard quantization of the classical PU oscillator 
based on the Hamiltonian (j4j) gives rise to a theory that is either non-unitary or unstable. They 
then take the dimcult-to-justify step of requiring that the coordinate z be considered as imaginary 
while the coordinate x, that is related to z via x — z, is treated as being real. This corresponds to a 
drastic change of boundary conditions that define the Hilbert space of the quantum theory. Using 
a similarity transformation they rotate the complex z-plane (z —> y := iz) to map the Hilbert space 
defined by this unusual boundary conditions to the usual Hilbert space L 2 (M. 2 ). However, in this 
Hilbert space the Hamiltonian takes the form 

which is manifestly non-Hermitian and PT-symmetric (the action of V being defined by x — > —x, 
Px — > —Pxi V y, and p y — > p y ). It turns out that (jHJ) is quasi- Hermitian [7], so that it can 
be mapped using another similarity transformation to a Hermitian Hamiltonian [8]. Using such a 
similarity transformation, one finds the Hamiltonian for an uncoupled pair of harmonic oscillators, 
namely H — 0^ + E^L^j _|_ + , that is clearly Hermitian and has a positive spectrum. 

The authors of [5] consider this procedure as a proof of a "No-Ghost Theorem for the Forth-Order 
Derivative PU Oscillator Model," though they admit that the quantum system they obtain does 
not have the classical PU oscillator (CQ) as its classical limit! 

The basic idea of the present investigation is that the forth-order equation of motion (JTJ that we 
use to define the classical PU oscillator may be derived using other Hamiltonians in the Hamiltonian 
formulation of classical mechanics. This suggests searching for a new classical Hamiltonian that 
generates the dynamical equation ([1]) and achieves the same goals as those of the approach of Ref. [5] , 
but avoids treating z as an imaginary variable. We will construct such a classical Hamiltonian for 
the non-degenerate PU Oscillator (where u\ 7^2) and elaborate on the degenerate case. 
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2 An Alternative Hamiltonian Formulation 



The starting point of our analysis is the introduction of a new variable, namely 

w ■= T 2 {z + Xz), (6) 

where r and A are a pair of nonzero free real parameters with dimension of time 2 and time -2 , 
respectively, so that w and z have the same dimension. It is an easy exercise to show that (CD) is 
equivalent to 

z = t~ 2 w — Xz, (7) 
w = (\- a )w-T 2 {X 2 -aX + (3)z. (8) 

We can view these equations as Newton's equations of motion for a system with two degrees of 
freedom. 

Next, we wish to obtain a Hamiltonian formulation of the dynamical system defined by ([7j) and 
(JSJ). As a first step we examine conditions on the free parameters r and A for which the system is 
subject to a conservative force. Introducing mass parameters /i 2 , fi w , and demanding the existence 
of a potential field V such that ([7]) and ([S]) take the form fi w iv = — d w V and fi z z = —d z V, we find 
by imposing the integrability condition, d w d z V = d z d w V, that 

A 2 - aX + (3 + = 0. (9) 

Introducing 

fi := (4/i-y ) 1/4 r-\ 5 := a 2 - 4/3 - fi 4 = (co 2 - u\f - fi 4 , (10) 
we can express the solutions of ([9]) in the form 

X = ^l. (11) 

With these choices for A, we can determine the potential field V, and introducing the momentum 
variables p w := fi w w and p z := fi z z, we obtain the following one-parameter family of classical 
Hamiltonians 

2 2 i 

Hs = ^ + ^ + ^(a T v^V 2 + ^(q ± V6)z 2 - -y/ '^ z (a 2 - 4/3 - S)wz. (12) 

Note that because we have not yet fixed the value of r and consequently Q, the parameter 5 is 
a free parameter. For all values of S, the Hamilton equations associated with Hs are equivalent 
to the the equation of motion ([T|) that defines the classical PU oscillator. A natural question is 
whether there are values of 5 for which the quantization of Hs would yield a stable and unitary 
quantum system. This turns out to require a separate analysis for the degenerate (u± = 0J2) and 
non-degenerate (u 1 ^ u 2 ) cases. 
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3 Non-degenerate PU Oscillator and Its Quantization 



As seen from (TT2i) . H$ is a real- valued function on the phase space provided that 5 > that is 
\f\oj\ — u;|| > f2. This is only possible for the non-degenerate case where u>i ^ u 2 . We will next 
show that in this case H$ is bounded from below. 

First we note that according to ([2]) and (ITOj) . a ± \/5 > 0. This allows us to introduce the real 
and positive parameters: 



and express f|T2|) as 



v w := \J jji w {a T VS), v z := -\hija ± v^), (13) 



rr . 2 2 , 2 2 2z/ w Z/ 2 fi 2 M72; 

pg, p 2 , / ^ 2 z V , 4/^V 



2/^ 2/x 2 V Va 2 - SJ a 2 -5 

In view of ([2]) and f lTUj) . the last term on the right-hand side of (1T4"1) is nonnegative. Therefore, H$ 
is bounded below by zero for all possible values of LO\ and uj 2 ^ toi (where 5 > 0). 

We can easily quantize the Hamiltonian ( 1T4"|) by applying the standard canonical quantization 
scheme: 

(w,z,p w ,p z ) (w,z,p w ,p z ), (15) 
where w, z,p w , and p z are operators acting in L 2 (R 2 ) as follows. For all ip G L 2 (IR 2 ), 

wip(w, z) = wip{w, z), zip(w, z) = zip(w, z), (16) 
p w ip(w, z) = -ihd w i/}(w, z), p z ij(w, z) = -ihd z tp(w, z). (17) 

The resulting quantum Hamiltonian operator, 

Hs - o h o h UwW / 9 f + "a 7' I 18 ) 

2/i M 2/i 2 V V « 2 - SJ ct ~ 

that also acts in the Hilbert space L 2 (M 2 ), is manifestly Hermitian, and being the sum of positive 
operators, it is a positive operator (with a non- negative spectrum.) It describes a stable and unitary 
quantum system consisting of a coupled pair of harmonic oscillators. By construction, taking the 
classical limit of this quantum system we recover the classical PU oscillator as defined by the 
equation of motion ([T]). 

Because (1T8]) is a quadratic Hamiltonian, we may try to decouple and diagonalize it using a 
linear canonical transformation. We defer the details to the appendix and suffice to mention that, 
for the non-degenerate PU oscillator that we consider here, there is a similarity transformation 
generated by a quadratic function Q of w, z,p w , and p z that maps H$ to the Hamiltonian operator 
H' s for a pair of decoupled harmonic oscillators: 

Hs^H>s = e^ 2 Hs e& 2 = A + f<^L + |L + i 2 , (19) 

where ii' w ,ji' z ,u± are positive real parameters whose explicit form are given in (1501) . below. 
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4 Degenerate PU Oscillator and Its Quantization 



For the degenerate PU oscillator, where 

U\ = uj<i =: (j, a = 2u 2 , 5 = — f2 4 , A = uj 2 ± — -, r = (20) 
the Hamiltonian ( Il2l reads 

H-n« = |^ + |^ + ^(2 W 2 T ^fi> 2 + ^(2 W 2 ±^fi 2 )^-iv^^^- (21) 
2(j, w 2/i z 4 4 2 

If we choose the mass parameters, that are actually not constrained by the classical equation of 
motion ([TJ, to coincide, i.e., set fi w = fi z =: /x, the Hamiltonian ( l2"Tj) takes the form 

#-Q 4 = ^(PL+Pl) + J pu J 2 T l ^)w 2 + {2uj 2 ± l tf)z 2 ]-^^. (22) 

Note that here Vt is a free nonzero real parameter. 

It is easy to see that regardless of the value of Q, the Hamiltonian H_q4 is invariant under the 
combined effect of complex-conjugation and swapping (w,p w ) and (z,p z ). Quantizing ([2~2l using 
the standard canonical quantization scheme ([TBI) , we obtain a non-Hermitian Hamiltonian operator 
(acting in L 2 (IR 2 )), namely 

= ^(p 2 +^ + ^[(2c, 2 T ^ (23) 

that possesses a particular antilinear symmetry. This is a generalized PT-symmetry whose gen- 
erator we denote by Ql~, [6]. Here Q is the operator of reflection about the line w = z in the 
w-z plane, and T is the time-reversal operator; for all ip G L 2 (M 2 ), (Qt/j)(w,z) = i/j(z,w) and 
(Tip)(w,z) = i/j(w,z)*. 

A major difference between our Hamiltonian formulation of the non-degenerate and degenerate 
PU oscillators is that in the latter case the classical Hamiltonian (121 j) is a complex-valued function 
of the phase space variables (w, z,p w ,p z ). Because Q is not constrained, we can take fiCu and try 
to use perturbation theory to see if the the spectrum of the Hamiltonian is real, j9]. Again, 

the fact that H-n A is a quadratic this Hamiltonian suggests the possibility of its diagonalization via 
a (possibly complex) linear canonical transformation. As we show in the appendix this turns out 
not to be possible. Indeed the lack of such a canonical transformation might be an indication that, 
similarly to the coupled oscillators studied in [U], is actually non-diagonalizable. We shall 

not pursue the study of the spectral properties of this Hamiltonian operator here, because as we 
explain below it does not define a quantum system that has the classical degenerate PU oscillator 
as its classical limit. 

Note that in view of (|2"U|) and (|2l), the parameter A and the coordinate variable w actually take 
complex values. This observation, which is often missed or ignored in the study of the classical 
systems underlying PT-symmetric quantum systems [UJ , reveals a basic deficiency of a quantiza- 
tion scheme involving w — > w, simply because while w is a complex variable the operator w, that 
is defined in f TT5|) . has a real spectrum. Performing the quantization of (12~T|) using ffTo^l does yield 
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a non-Hermitian Hamiltonian operator H_qa with an antilinear symmetry, but even if this Hamil- 
tonian turns out to be quasi-Hermitian and capable of defining a unitary quantum system, this 
system does not admit the classical PU oscillator as its classical limit. The same conclusion applies 
to the analysis of [5]. The main difference is that in our approach this problem only arises for the 
degenerate PU oscillator, while in the approach of Ref. [5] it is present also for the non-degenerate 
PU oscillator. 

Next, we recall that in view of O20p . Eqs. (JT]) , (EJ) , and OH]) take the form 

z< 4 > + 2u 2 z^ + u 4 z = 7 (24) 
/iz = kw — (k ± i&)z, (25) 
fiw = (—k ± i&)w + kz, (26) 

where k := /icu 2 and ^ := ^y-. Note that, as in the non-degenerate case, here z is a real-valued 
function of time, whereas w is necessarily complex-valued. This is because (l25|) and ( |26l) yield ( |24l 
only for 8. ^ 0. 

We can express Eqs. (1251) and ( 1261) in terms of real and imaginary parts of w, namely w\ := 
Re(u') and w 2 := Im(w). These yield a system of three Newton's equations for z, w\ and w 2 that 
can be reduced to 

fiz = — /c2, = —kwi, w 2 = ±z. (27) 

It is not difficult to show that the force corresponding to the first two of these equations is non- 
conservative. This means that this system does not admit a Hamiltonian formulation based on a 
Hamiltonian that is quadratic in momenta. In other words, we can obtain a consistent Hamilto- 
nian formulation of the classical dynamics of the system using the Hamiltonian ff22l . only if we 
treat both the coordinate variables w and z and the corresponding momentum variables p w and p z 
as complex variables. This yields an analytic continuation of the degenerate PU oscillator into a 
complex dynamical system defined by ( 1221) via the complex Hamilton equations 

w = d Pw H, z = d p M, p z = -d z H, p w = -d w H, (28) 

where we use H in place of H_^a for simplicity. 

A detailed study of complex classical systems of the type (1281) (though with a single complex 
coordinate variable) has been carried out in [12]. Equations (1281) are meaningful if H is a complex 
analytic function of w, z, p w , and p z . This implies that Cauchy-Riemann conditions hold for the 
real and imaginary parts of H. That is, introducing 

zi := Re(z), z 2 := lm(z), p Zl := Re(p z ), 

p Z2 := Im(p 2 ), p Wl := Re(w z ), p W2 := lm(w z ), 

H 1 := Re(if), H 2 := lm(H), 

d Wl Hi = d W2 H 2 , d W2 H\ = —d Wl H 2 , (29) 

d Pwx H x = d Pw2 H 2 , d Vw H x = -d Pwi H 2 , (30) 

8 Z1 H X = d Z2 H 2 , d Z2 H x = -d Zl H 2 , (31) 

d Pzi H 1 = d Pz2 H 2 , d Pz H x = -d Pz E 2 . (32) 



we have 
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As originally noted in [12] for one- dimensional complex configuration spaces, the following are 
miraculous consequences of the Cauchy-Riemann conditions (12"5|) - (1521) . 

• Eqs. (|28|) are equivalent to a system of Hamilton's equations, ij = d Pi M ', pj = —d Xj J^, for 
the phase space variables 

x 1 :=w 1 , x 2 =p W2 , x 3 :=z 1 , x A := p Z2 , 

Pl-=P W1 , P2 = W 2 , P3-=P Z1 , Pi-=Z 2 , 

where M' is nothing but the real part of the complex Hamiltonian H, i.e., Jif := H±. 

• The imaginary part J? := H 2 of the complex Hamiltonian H is an integral of motion for the 
system that is independent of J^ 3 , 

4 

{^,^} := [(d Xj y)(d Pj Jff) - (d Pj J?)(d xj Jf)] = 0. 

3=1 

The problem of quantizing one-dimensional examples of similar complex dynamical systems 
has been considered in [13]. The integral of motion J? generates a type of gauge symmetry for 
these systems [T21 [T3] that makes their quantization a nontrivial and non-unique process. Here 
we will only derive the explicit form of and J? to examine if a canonical quantization of J^ 7 
(before imposing the constraint) can give rise to a stable and unitary dynamics. A straightforward 
calculation gives 

& = ^^-^+^)+^4±^(zi^-x 3 P4)-^^ + |(^ + ^)-^a:iX3 J (33) 

= i-g (P 2 - Pa) H 1" «(ZlP2 + ^3P4j - .£(ZlP4 + ^3P2j ± ^( X 3 ~ I 34 ) 

As seen from (1331) . ^ is an unbounded function of the phase space variables. Therefore, the quan- 
tization of J$? leads to a Hermitian operator jj? whose spectrum is not bounded below. However, 
note that to impose the constraint $ := J' — C = 0, where C is a real constant, we need to quantize 
<3> to obtain the quantum constraint $ and define the physical Hilbert space of the system as the 
null space (kernel) of $. Given the complicated form of J^f and it is not easy to determine if 
the dynamics taking place in the physical Hilbert space is stable. But this seems to be unlikely. 

Alternatively we can impose the constraint before quantization by solving $ = for one of the 
variables and substituting the result in the expression for Jif. This yields a highly complicated 
reduced classical Hamiltonian that is still unbounded both from above and below. Hence its 
quantization gives rise to an unstable quantum system. 

5 Concluding Remarks 

We have established a novel Hamiltonian formulation of the Pais-Uhlenbeck oscillator that yields 
a real and positive classical Hamiltonian for a non-degenerate oscillator. Canonical quantization of 
this Hamiltonian yields a unitary and stable quantum system with the classical PU oscillator as its 
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classical limit. This provides a simple and consistent solution of the long-standing non-unitarity 
versus instability problem that one encounters in quantizing the classical PU Oscillator. 

Our Hamiltonian formulation involves a manifestly complex classical Hamiltonian whenever the 
PU oscillator is degenerate. The naive canonical quantization of this Hamiltonian, that allows 
for mapping complex classical variables to Hermitian operators, gives rise to a non-Hermitian 
Hamiltonian operator that possesses a generalized PT-symmetry. We are of the opinion that even 
through this Hamiltonian has appealing symmetry properties, it is not relevant to the problem 
of quantizing the classical PU oscillator. This is because this approach is void of a consistent 
quantum-to-classical correspondence that would map classical observables to quantum observables. 

A consistent method of dealing with this complex classical Hamiltonian is to realize that it ad- 
mits a real description. Using this description one discovers that the dynamics is actually generated 
by the real part of the complex Hamiltonian H and that the system has a highly nontrivial 
integral of motion given by the imaginary part J" of H. The latter may be identified both as the 
generator of a gauge symmetry and as a first class constraint. A proper quantization of this system 
must take into account the presence of this constraint. The Hamiltonian operator 34? that one 
obtains by quantizing 3f is a Hermitian operator with an unbounded spectrum both from above 
and below. In principle this is not sufficient to conclude that the quantum dynamics is unstable, 
because the physical Hilbert space "H p h ys is the null space of the quantum constraint $ = J> — C. 
To determine the stability of the quantum system after the imposition of the constraint, one must 
study the spectrum of the restriction of Jt? to "H p h y s- Given the complicated form of 3if and $ this 
is not an easy task. There are however hints that the degenerate quantum PU oscillator probably 
remains unstable even after enforcing the constraint. The presence of runaway classical solutions 
seems to support this conclusion. 

Acknowledgment: This work has been supported by the Turkish Academy of Sciences (TUBA). 



Appendix 

In this appendix we explore the possibility of decoupling and diagonalizing the Hamiltonian operator 
obtained by standard canonical quantization (|15p of (I12p. namely 

H& = ^ + ^ + ^f{a =F V6)w 2 + ^(a± V6)z 2 - \y/ ^{a* - 4/3 - 5)wz. (35) 
First, we introduce the dimensionless quantities: 



h 



w, x 2 :-- 



M/« ; 

h 



._,V6_, VW-uiy-w 

^ " 2 i 2 ' 



-z, p! := 



7 - 



P2 



Pz 



a 



(36) 
(37) 



that allow us to express Hg in the form 



X\X 2 



(38) 
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It is easy to check that = iSij for i,j = 1,2, and a/7 2 — e 2 = fi 2 /a > 0. 

A linear canonical transformation is equivalent to the similarity transformation Hs — >• H' := 
e -Q/ 2 ff (5e Q/ 2 where Q is a quadratic function of the operators Xj and with z = 1, 2. It is well- 
known that the transformed operators Xj := e~®/ 2 Xie®/ 2 and Pj := e~^l 2 pie^l 2 are linear functions 
of £j and p, satisfying the canonical commutation relations 



(Pi,x 1 ,p2,x 2 ) T and q' := 



iSij. (39) 

Let us introduce q := (pi, x±,P2, Z2)' 1 and q' := (Pi,Xi, P 2 ,X 2 ) T , where the superscript "T" 
stands for the transpose. Then we can express and H' s in the form 

q T H 5 q, 



q /T H 5 q', 



where 



H 



8 - z 



( 1 













1 












\ 



J 



(40) 



(41) 



The linearity of the canonical transformation is equivalent to the existence of an invertible 4x4 
complex matrix U (with numerical entries) such that 



q = Uq. 

In terms of U the canonical commutation relations ( 1391) take the form 

U T CU = c, 

where 



/0 
1 




1 












1 






-1 




(42) 



(43) 



(44) 



We wish to find out if we can choose Q such that H' s is the Hamiltonian operator for a pair of 
decoupled oscillators. In view of ( 140]) . this is equivalent to looking for a matrix U such that the 
matrix 

:= U T H 5 U (45) 

is diagonal. This problem has been studied in [15]. Here we outline its solution. 
Using ( T4"5l) and the properties of C, we find that 



(CH 5 ) 2 U = u(-ch^c t h^; 



(46) 



Now, for a diagonal H^, the matrix — CH' 5 C T H^ is also diagonal, and (T46j) implies that the columns 
of U are eigenvectors of (CH^) 2 . We can easily compute and solve the eigenvalue problem for this 
matrix and construct U. The result is 

/ 



U 






e+7 





e-7 


\ 






<E+7 





e-7 







V-7 2 -* 2 









1 





1 




1 





1 





/ 



(47) 
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Substituting (I47p in (145]) . we have checked that indeed is a diagonal matrix with diagonal 
entries: 



hy/a^jl - 7) 
2(7 - e) 



(5J22 



hy/aj 



(Hi 



7 



V«7(l+7) 
2(7 + e) 



(Hi) 



ii 



h\fa 7 
7 + e 



• (48) 



This implies that H' s describes the dynamics of a pair of decoupled oscillators: 



q' T lW = q T Hiq 



hy/a 



7(1 -l)Pj 2 7 fj 7(1 + 7)^2 



7 



7 - e 



7 



2tx| 
7 + e 



where 



Pw 

2/4 



- U> 



(7 - 



24 



+ /M,2 



// 2 



(7 + e)^ 



V2 7 2 (l±7)« 



(49) 



(50) 



7(1-7)' 7(1 + 7)' 7±e 

For the non-degenerate PU oscillator where 5 is a nonnegative real number, fi' w , fi' z and u;± are 
real and positive parameters, and H' 6 is the sum of two noninteracting simple harmonic oscillators. 
Clearly the spectrum is given by {/i[u;_(n_ + §) + u + {n+ + -)] | n± = 0, 1, 2, • • • }. This coincides 
with the spectrum of i^, for if 5 and H' s are isospectralJll 
Next, we note that according to ( 14"7|) . 



detU 



4 7 2 



(51) 



Therefore, for a degenerate PU oscillator where 7 = 0, there is no invertible matrix U that could 
implement the desired linear canonical transformation. This means that in this case we cannot 
decouple and diagonalize the Hamiltonian H$ by performing a linear canonical transformation (a 
similarity transformation with a quadratic generator Q.) This suggests that in this case Hg might 
even be non-diagonalizable|f| 



1 In fact, in this case the similarity transformation can be implemented using an anti-Hcrmitian generator Q. 
Therefore it is a unitary transformation that leaves the spectrum invariant. 
2 We are currently unable to establish if this is actually true. 
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